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Within the framework of Dyson-Schwinger equations (DSEs), we discuss the chiral phase transi¬ 
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of a scalar susceptibility, we find that there might be no CEP 5 in the T — /is plane, and the phase 
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Key-words: chiral chemical potential, chiral phase transition, Dyson-Schwinger equations 
PACS Number(s): ll.30.Rd, 25.75.Nq, 12.38.Mh, 12.39.-x 


I. INTRODUCTION 

Quantum Chromodynamics (QCD), which describes 
the interactions between quarks and gluons, is already 
commonly accepted as an essential part of the Standard 
Model of particle physics. Dynamical chiral symmetry 
breaking (DCSB) and quark color confinement are two 
fundamental features of QCD, and there are also many 
laboratories and experiments on this field, such as the fa¬ 
mous Relativistic Heavy Ion Collider (RHIC) and Large 
Hadron Collider (LHC). However, thanks to the com¬ 
plicated non-Abelian feature of QCD itself, it is so dif¬ 
ficult to have a thorough understanding of the mecha¬ 
nisms of DCSB and confinement, especially in the in¬ 
teresting non-perturbative region, which means quarks 
and gluons are strongly coupled to each other and then 
the related processes have small momentum transfer (or 
equivalently, the coupling constant becomes large and 
running). In this case, nowadays people often and in 
some sense have to resort to various effective models to 
study them phenomenologically, such as the chiral per¬ 
turbation theory [1-5], the global color symmetry model 
(GCM) [6-10], the quasiparticle model [11-18], the QCD 
sum rules [19-22], the Nambu-Jona-Lasinio (NJL) model 
and the related Polyakov-loop-extended Nambu-Jona- 
Lasinio (PNJL) model [23-34], Lattice QCD [35-37], and 
the Dyson-Schwinger equations (DSEs) [38-45] . Through 
these studies, people hope to get profound insight of our 
nature as well as the early Universe. 


* Email: zonghs@nju.edu.cn 


Generally speaking, chiral symmetry is an exact global 
symmetry only when the current quark mass m is zero 
(the chiral limit case). In the low temperature (T) and 
low chemical potential (/i) phase (the hadronic phase, 
often referred to as Nambu-Goldstone phase or Nambu 
phase), this symmetry is spontaneously broken, and as 
a consequence there exist — 1 {Nf is the number of 
flavour) pseudoscalar Nambu-Goldstone bosons, mean¬ 
while the QCD vacuum hosts a chiral condensate (two 
quark condensate) (qq) (which can actually act as an or¬ 
der parameter for chiral phase transition). At present, 
it is commonly accepted that when temperature and/or 
quark chemical potential are high enough, the strongly 
interacting hadronic matter will undergo a phase tran¬ 
sition to some new phase, where the chiral symmetry 
is restored for the chiral limit case or partially restored 
for the m ^ 0 case. This new phase is usually called 
Wigner phase, and in some sense is related to the fa¬ 
mous quark gluon plasma (QGP), which is expected to 
appear in the ultra-relativistic heavy ion collisions or the 
inner core of compact stars. As for the nature of the 
chiral phase transition when to 7^ 0, a popular scenario 
favors a crossover at small chemical potential, and then 
turning into a first order chiral transition for larger chem¬ 
ical potential at a critical end point (CEP) [46]. Such a 
picture is consistent with most Lattice QCD simulations 
and various QCD-inspired models, as listed in the last 
paragraph, however, it is not yet clarified directly from 
the first principles of QCD. The search for such a CEP is 
also one of the main goals in the high energy physics ex¬ 
periments, such as the beam energy scan (BES) program 
at RHIC [47-50]. Unfortunately, Lattice Monte Carlo 
simulations cannot be used to resolve this issue directly 



2 


due to the “sign problem” [51-56], and until now there is 
still no firm theoretical evidence for the existence of such 
a CEP, so the calculations based on some effective QCD 
models are also irreplaceable nowadays. 

In Ref. [57], K. Fukushima et al firstly introduce the 
chiral chemical potential (also called axial chemical po¬ 
tential in some other literatures, such as Ref. [58]), /is, 
which is conjugated to chiral charge density; and in 
Ref. [59], M. Ruggieri suggests that the CEP of the chiral 
phase diagram can be detected by means of Lattice QCD 
simulations of grand-canonical ensembles with this chiral 
chemical potential. By concrete calculations within some 
chiral models, the author shows that a continuation of the 
CEP at finite temperature and finite chemical potential, 
to a possible “CEP 5 ” in the T — /rs plane is reachable, 
which is then helpful in the determining of the CEP in 
the T — ^ plane from Lattice QCD. The existence of such 
a possible CEP 5 is also confirmed in some other chiral 
model studies, for example, Refs. [60, 61]. In Ref. [58], 
the authors investigated the effect of the vector inter¬ 
action as well as the finite current quark mass on the 
location of the CEP. In this paper, we will discuss the 
related topics within the framework of Dyson-Schwinger 
equations, which is widely used as well as has been proved 
to be successful in hadron physics and phase transitions 
of strongly interacting matters. The following of this pa¬ 
per is organized in such a way: in Sect. II we give a 
basic introduction to the DSEs at finite temperature and 
nonzero chemical potential as well as an effective model 
gluon propagator, and with the help of a scalar suscepti¬ 
bility we also discuss the nature of the chiral phase tran¬ 
sition within this framework; then in Sect. Ill, we discuss 
the influences of the chiral chemical potential on the chi¬ 
ral phase transition of QCD in detail, and mostly focus 
on the behaviours of the CEP, not only the algebra but 
also the numerical results; at last, a brief summary is 
given in Sect. IV. 


II. DYSON SCHWINGER EQUATIONS AND 
AN EFFECTIVE GLUON PROPAGATOR 

In this section, we will briefly review the formula of 
Dyson Schwinger equations, which is not only widely 
used in the non-perturbative region of QCD, but also 
in some other fields like the Quantum Electrodynamics 
in (2-1-1) dimensions (QED 3 ) [38, 62-65], etc. At zero 
temperature and zero chemical potential, the DSE for 
the quark propagator reads [38] (we will always work in 
Euclidean space and take the number of flavors Nf = 2 
while the number of colors Nc = 3 throughout this pa¬ 
per. Moreover, as we employ a ultra-violet finite model, 
renormalization is actually unnecessary) 

S{p)-^ = So{p)~^ + ^J 

( 1 ) 


where S{p) is the dressed quark propagator, 

So{p)~^ = h ■ p + m, (2) 

is the inverse of the free quark propagator, g is the strong 
coupling constant, D^^{p — q) is the dressed gluon prop¬ 
agator, and Pjy = T^{p, q) is the dressed quark-gluon ver¬ 
tex. According to the Lorentz structure analysis, we have 

5(p)-i=*Mp") + R(/), (3) 

where A{p^) and B{p^) are scalar functions of p^. Af¬ 
ter the gluon propagator together with the quark-gluon 
vertex are specified, people can then solve this equation 
numerically. 

The extension of the above quark DSE to its nonzero 
temperature and nonzero quark chemical potential ver¬ 
sion is systematically accomplished by transcription of 
the quark four-momentum via p Pk = where 

Wfc = ojk + ip with Wfc = {2k + k G Z the fermion 

Matsubara frequencies, and no new parameters are in¬ 
troduced [39] 

Sipk)~^ = So{pk)~^ + ^T'^g'^Df.^ipk - qnhk.^i<ln)T^<^- 

(4) 

where 

So{Pk)~^ = ij ■ Pk + m, ( 5 ) 

and^ denotes I ■ Nevertheless, its solution 

now should have four independent amplitudes due to the 
breaking of 0(4) symmetry down to 0(3) symmetry [39] 

S(pk)~^ =i fA{p‘^,[bl) + lR(p^w^) 

-f^74a;feO(p^w^)-k f-fiCokD{p'^ (6) 

where ^ = 7 • p, 7 = ( 71 , 72 , 73 ), and the four scalar 
functions T = A, B, C, D are complex and satisfy 

B{p\Qlr =B{p\i:jl,_,). (7) 

But as discussed in Ref. [39], the dressing function D is 
power-law suppressed in the ultra-violate region, so that 
actually does not contribute in all cases investigated in 
our work. At zero temperature but nonzero chemical po¬ 
tential case, D vanishes exactly since the corresponding 
tensor structure has the wrong transformation proper¬ 
ties under time reversal [66]. For these reasons, in most 
cases we can just neglect D, and get the commonly used 
general structure of the inverse of quark propagator as 

S{pk)~^ = i /A(p^w^)-hlR(p^w^,)-h^74WfeC'(p^w^). 

( 8 ) 

For the dressed-gluon propagator, the general form is 
like this, 

g^D^^{kni) = P]^^Dt{P + pI;^Dl{P (9) 

where kni = {k,ujni) = (p- g,a;„ - wj), are trans¬ 
verse and longitudinal projection operators, respectively. 
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And for the domain T < 0.2 GeV, of which we are con¬ 
cerned in this work, the authors of Ref. [67] have proved 
that Dt = Dl is a good approximation. For the in¬ 
vacuum interaction, in this work we will adopt the fol¬ 
lowing form of Ansatz as in Ref. [68], 

= = ( 10 ) 

(T° 

which is a simplified version of the famous as well as 
widely used one in Refs. [69, 70]. It can be proved that 
this dressed gluon propagator at T = 0 violates the ax¬ 
iom of reflection positivity [71], and is therefore not ob¬ 
servable; i.e., the excitation it describes is confined. The 
same is true of the dressed quark propagator which is also 
not positive definite and hence is confined (Actually, we 
can take the gluon propagator as input, and the quark 
propagator can then be solved numerically. The results 
show that there is no singularity on the real, positive, i.e. 
timelike, axis, which implies that quarks are confined). 

As concerning the quark-gluon vertex, in this work we 
will take the rainbow truncation, which means a simple 
but symmetry-preserving bare vertex is adopted, 

= ( 11 ) 


The status of propagator and vertex studies can be 
tracked from Ref. [72]. 

Now let us fix the related parameters and then show 
some of the numerical results. Dq and a are usually 
fixed by fitting the observables, such as the two-quark 
condensate, the pion decay constant (/tt = 131 MeV) 
and the pion mass (m.^. = 138 MeV). In this work we 
adopt the one from Ref. [70], that Hq = 9.3 x 10® MeV^ 
and cr = 400 MeV. For the current quark mass we will 
use to= 5 MeV. Then substituting Eqs. (5), (8), (10), and 
(11) into Eq. (4), we can solve the quark DSE for each 
value of temperature and chemical potential by means of 
numerical iteration. As an example, we show 5(0, Wg) as 
a function of /i for different T in Fig. 1, and the corre¬ 
sponding chiral susceptibility with respect to m, which is 
defined as 


Xm(T, p) 


dm 


( 12 ) 


in Fig. 2. 

In general, we can see from Fig. 1 that the scalar func¬ 
tion 5(0, ujq) will decrease when the chemical potential /r 
increase, this phenomenon holds to be true for the tem¬ 
perature T and momentum too. It is known that the 
scalar part B{p‘^,Col) of the quark propagator Eq. (4) in 
some sense reflects the dressing effect of the quark, so the 
results show that the dressing effect becomes weaker and 
weaker for higher T, p and . We can also see from Fig. 1 
that for different values of T, 5(0, Wg) may behave differ¬ 
ent: for T larger than a critical Tc = 129 MeV, 5(0, Wg) 
change gradually but continuously from the Nambu solu¬ 
tion to the Wigner solution; while for T smaller than T^, 
there will appear a sudden discontinuity at some critical 

p. 
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FIG. 1. B(0,luo) as a function of p for three different T. 
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FIG. 2. Xm(r, p) as a function of p for three different T. 


To study the nature of the chiral phase transition, es¬ 
pecially to determine the critical value of pc at Tc, peo¬ 
ple often employ various susceptibilities of QCD [73-77]. 
We can see from Fig. 2 that, for T > Tc, the suscepti¬ 
bility Xm indicate a crossover from the Nambu phase to 
the Wigner phase, and the peak grow higher and higher 
when T approaches Tc- At Tc, Xm shows a sharp and 
narrow divergent peak, and the value of this peak turns 
to be oo, which demonstrate that here is a second-order 
phase transition, and corresponding to the CEP. And for 
T < Tc, an obvious first order phase transition will occur. 
According to these results, we can move on to study the 
chiral phase transition, especially the behavior of CEP. 
In the following Sec. Ill, we will focus on the variance of 
CEP when the chiral chemical potential is considered as 
an additional scale. 


III. INFLUENCES OF THE CHIRAL 
CHEMICAL POTENTIAL ON THE CHIRAL 
PHASE TRANSITION 


The concept of chiral chemical potential was first pro¬ 
posed by K. Fukushima et al. in a study related to the 
external magnetic field [57]. Since topological charge 
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changing transitions can induce an asymmetry between 
the number of right- and left-handed quarks due to the 
axial anomaly, they introduce the chiral chemical po¬ 
tential /i 5 , which couples to the difference between the 
number of right- and left-handed fermions. The chirality 
is also expected to be produced in the high tempera¬ 
ture phase of QCD [58]. Many researchers argue that 
although /is is a mere mathematical artifice instead of a 
true chemical potential, it has the advantage that can 
be simulated on the Lattice QCD with N(, = 3, hence is 
likely to provide some useful information for the studies 
of the CEP, and even for inhomogeneous phases or the in¬ 
ner structure of compact stellar objects. One of the most 
interesting features of the introduction of /X 5 is that, it 
makes the continuation of the CEP to a possible “CEP 5 ” 
in the T — /is plane possible, which is of course helpful 
in the determining of the CEP from Lattice QCD. Some 
other researchers also confirm such a possible CEP 5 in 
related chiral model studies [60, 61]. In this part, we 
will discuss the related topics within the framework of 
Dyson-Schwinger equations. 

To be specific, in order to study the effects of /is, people 
should add the following term to the Lagrangian density, 

M5'07475'0- (13) 


Por the details of S')®, /is), please see the Appendix part. 

Substituting Eqs. (15), (9), and (11) into Eq. (17), we 
found that the solution is the following coupled integral 
equations (for the sake of concise, here all the notion J^ip) 


means T(di 




A(p) 


khe-^- 

X Ka. 

B(p) 

= m-kc(r)^ 

khe-^ 

x/Cs, 

C{p) 

= l + c(T)^, 

^2 -kl 

X /Cc, 

D(p) 



X ICd, 

M{p) 

-ciTp.e 


X /C 4 I 5 , 

B^{p) 



X /Cbs, 

d(p) 

= f/xs/w„c(r)] 

^khe 

. 2/2 

x/Cc, 

Dz{p) 

= c{T)'tkl,e 


X /Cdj, 


in which 


And in our work, the quark propagator and its inverse 
now can include at most eight components according to 
the Lorentz structure analysis, namely, 

1, f, 74, #74, 75, #75,7475, #7475- (14) 

So now the general inverse form of the dressed quark 
propagator is 

S{Pn, = i fA + B + i 74 W„C'-|- ^ 4 ,u}„D 
-|-(i /As + B 5 + ijiu;nC5+ /74<l’nDs)7s, (15) 

The eight scalar functions B = A, B, C, D, As, i?s, Cs, 
Ds denotes /is), which are all complex and 

satisfy the following equation, 

= J^(/^W?.fe_l,/is)• (16) 

Now the quark DSE at nonzero temperature and 
nonzero chemical potential is then, 

S{Pn,Pb)~^ = So{pn,Pb)~^ 

- m)ltiS{quPf,)Tv{pn,qi), (17) 

where 

So{Pn,P5)~^ = i f+m + f74W„ - /is7475- (18) 


ICa = -[(/■ qkni + 2k-pk ■ ^ua + 2k ■ pLOniuJiac]/i'j^kli), 

JCb = 5a b, 

JCc = -uJiaduJn - 2uJni{k ■ qaA + UJnlUJiac)/{uJnkli), 

K^d = [p- - kh) + 2k ■ pk 

= [{p-qkh +2k ■ pk ■ ^aAp, + 2k ■ p^nmcrd/il^kh), 
5aB^j 

= UJiac/uJn+2uJni{k ■ qo-Ap, +U!nlUJiaCp,)/iu!nkli), 

= [p-qi‘^‘^li-kli)+2k-pkd>^i<^Dj{fi^nkli). (20) 

Then, we can solve Eq. (19) numerically for specific chi¬ 
ral chemical potential /is, as well as for the temperature 
T and normal chemical potential /i. The critical Tc[ph] 
and /ic[/is], which are coordinates of the new “critical 
end point” at a specific /is (CEP [/is]), are determined 
by the scalar susceptibility that is defined in Eq. (12) 

In this work, we will concentrate on the behavior of the 
CEP [/is], which is expected to be linked to a possible 
CEPs in the T — /is plane [59]. In Pig. 3 and Pig. 4, we 
plot the relations between /is and 7), [/is] as well as /is 
and /ic[/is], respectively, which are obtained by seeking 
the corresponding CEP [/is] for different /is. Therefore, 


The reason for this is easy to understand, since the difference in 
densities of the right- and left-handed quarks, 77,5 = nji — ul, 
does not conserve. 


Here we summarize that, in this work “CEP” means the critical 
end point in the T — ^ plane with /is = 0, while CEP 5 denotes 
the possible one in the T —/is plane with /i = 0 [59], and CEP [/is] 
is generally the similar critical end point in the T — plane for 
a specific /is- 
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FIG. 3. The relation between /rs and the corresponding Tcl/rs] 
in the T — ^ plane. 



Ml^s] (MeV) 


FIG. 4. The relation between /rs and the corresponding 
in the T — ^ plane. 


each point in these two lines means a “CEP’ in the T — ^ 
plane for the corresponding 

We can see Fig. 3 that, Tc[n 5 ] increase slowly when 
fj ,5 is smaller than approximately 300 MeV, and turns to 
increase quickly for larger fj,^. Nevertheless, the whole 
increase is smooth and monotonic. The most interest¬ 
ing thing is as shown in Fig. 4, that /tci/ts] will decrease 
firstly when ^5 is smaller than about 400 MeV, and then 
increase for larger The behavior of Hc[fJ- 5 ] for differ¬ 
ent ^5 is quite different with the previous results from 
some chiral models, such as Fig. 4 of Ref. [59], that at 
some critical value of /xs, fJ,c[fJ- 5 ] will decrease to 0 , where 
a CEP 5 is expected to exist. When ^5 is not very large, 
the qualitative properties of the results in Fig. 4 are simi¬ 
lar to those from chiral models, but the decrease of ttciMs] 
is much slower. Our results also indicate there might be 
no CEP 5 in the T —^5 plane, and the phase transition in 
the r — ^5 plane might be totally crossover when ^ < 50 
MeV, which then has apparent consistency with the Lat¬ 
tice QCD calculation [53]. Furthermore, our studies may 
also provide some hints for recent studies related to the 
chiral chemical potential, such as Ref. [78]. Since com¬ 
paring with the chiral models, DSEs are renormalizable. 


and include the effects of color confinement as well as 
DCSB simultaneously, accordingly the DSEs approach 
is commonly accepted to be a closer theory to QCD it¬ 
self [38-41, 43]. In the calculations within the framework 
of DSEs, people do not need to introduce some annoying 
parameters, such as the momentum cutoff scale which 
destroys some basic symmetries of QCD, and then can 
give more reliable results. In other words, although the 
chiral models have been found to give reasonable phase 
diagrams hiT — ji plane, they can not guarantee the ro¬ 
bustness when jjL^ acts as an additional scale, and the 
introducing of /is might make the applicability of the 
calculations T -b /i -I- /is < A (A is some cutoff scale). 

IV. SUMMARY 

Thanks to the complicated non-Abelian feature of 
Quantum Chromodynamics (QCD) itself, its two fun¬ 
damental features, namely, dynamical chiral symmetry 
breaking (DCSB) and quark color confinement, have to 
be studied phenomenologically through various effective 
models at present, especially in the most interesting non- 
perturbative region. In this work, we discuss the chiral 
phase transition of QCD within the framework of Dyson- 
Schwinger equations (DSEs), with a chiral chemical po¬ 
tential /is as an additional scale other than the normal 
temperature T and quark chemical potential /r. We give 
a basic introduction to the DSEs at finite temperature 
and nonzero chemical potential as well as an effective 
model gluon propagator firstly, and then mainly focus on 
the calculations related to the famous critical end point 
(CEP) in the T — /r plane, which is predicted by many 
model studies, and has caused much interests both in the 
experimental side (one of the main goals in some high 
energy physics experiments) and theocratical side. With 
the help of a scalar susceptibility, that often act as an 
order parameter of chiral phase transition, we find that 
there might be no CEP 5 in the T — /rs plane, which is 
thought to exist by some chiral model calculations, and 
the phase transition in the T — /xs plane might be to¬ 
tally crossover when /r < 50 MeV, which has apparent 
consistency with the Lattice QCD calculation. DSEs is 
widely used as well as has been proved to be successful 
in hadron physics and phase transitions of strongly inter¬ 
acting matters, so that our study may also provide some 
useful hints to some other studies related to //s. Last but 
not least, we’d like to say that the related issues deserve 
further studies. 
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APPENDIX: STRUCTURE OF THE QUARK 
PROPAGATOR 

Using Eq. (15), after some algebra we find that the 
quark propagator can be written as (note that here T = 

S{qhiJi.b)=i icy A + (yB + i^iCbi(Jc+ |74Wi ctd 

+(* i<yA^ + CJB^ + Z74Wicrc6+ ijiUiiao^h^, (21) 


where 

CA = { — Ati + C^UIlt2) 

(JB = {Bti - D^ujift2)/h, 

oc = {-Cti - A^^t 2 lu!i)IU, 

CCD = {—Dti — B5t2/u;i)/t3, 

(JA^ = {A^ti - CuJlt2)/t3, 

UB^ = {-B^ti + Du}ii^t2)/t3, 

O-Cs = (-C'sU - t2M)/t3, 

— Bt2juJi)jt3, ( 22 ) 


ti=B^-Bl + (C^ - Ci)iof 

+ [A^.Al + {D^-DlU]cf, 
t 2 = 2wi{A5C + B 3 D - AC 3 - BD 3 ), 
t3 = tl+cftl, (23) 
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